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ABSTRACT 

We present a consistent three-dimensional model for the head-tail radio galaxy NGC 1265 that 
explains the complex radio morphology and spectrum by a past passage of the galaxy and radio bubble 
through a shock wave. Using analytical solutions to the full Riemann problem and hydrodynamical 
simulations, we study how this passage transformed the plasma bubble into a toroidal vortex ring. 
Adiabatic compression of the aged electron population causes it to be energized and to emit low 
surface brightness and steep-spectrum radio emission. The large infall velocity of NGC 1265 — which 
is barely gravitationally bound to the Perseus cluster at its current position — and the low Faraday 
rotation measure values and variance of the jet strongly argue that this transformation was due to 
the accretion shock onto Perseus situated roughly at i?2oo- Estimating the volume change of the 
radio bubble enables inferring a shock Mach number of ~ 4.2t?;2i ^ density jump of 3.41q;4, a 
temperature jump of 6.3^2 ^^id a pressure jump of 21.5 ± 10.5 while allowing for uncertainties in the 
equation of state of the radio plasma and volume of the torus. Extrapolating X-ray profiles, we obtain 
upper limits on the gas temperature and density in the infalling warm-hot intergalactic medium of 
kT < 0.4 keV and n < 5 x 10~^cm~'^. The orientation of the ellipsoidally shaped radio torus in 
combination with the direction of the galaxy's head and tail in the plane of the sky is impossible to 
reconcile with projection effects. Instead, this argues for post-shock shear flows that have been caused 
by curvature in the shock surface with a characteristic radius of 850 kpc. The energy density of the 
shear flow corresponds to a turbulent-to-thermal energy density of 14% — consistent with cosmological 
simulations. The shock-injected vorticity might be important in generating and amplifying magnetic 
fields in galaxy clusters. We suggest that future polarized radio observations by, e.g., LOFAR of head- 
tail galaxies can be complementary probes of accretion shocks onto galaxy clusters and are unique in 
determining their detailed flow properties. 

Subject headings: galaxies: clusters: individual (Perseus) - galaxies: individual (NGC 1265) - galaxies: 
jets - intergalactic medium - radio continuum: galaxies - shock waves 



1. INTRODUCTION 

Head-tail radio galaxies show spectacular asymmetric 
radio morphologies and occur in clusters of galaxies. The 
favored interpretation of these head-tail sources is radio 
jets ejected by an active core of a galaxy. The jets are 
bent at some angle toward one direction giving rise to 
a "head" structure and fan out at larger distances in a 
characteristic tail that extends over many tens to hun- 
dreds of kpc. If the flow impacting these galaxies is su- 
personic, ram pressure of the intracluster med ium (ICM) 
causes the jets to bend (jBegelman et al.l[l979t ): if the flow 
is trans-sonic with a Mach number 7V1 ~ 1, the ther- 
mal pressure gradient of the interstellar medium of these 
galaxies due to their motion t hrough the ICM causes the 
bending (jJones fc OwenlflQT^) . In this model, the super- 
sonic inflow is decelerated and heated by a bow shock in 
front of the galaxy which also generates a turbulent wake 
that re-accelerates the relativistic particle population in 
the tail and illuminates the large-scale tail structure of 
these sources. The Perseus cluster is an outstanding ex- 
ample, as it hosts many of these head- tail sources: most 
prominently NGC 1265 and IC 310 (jRvle fc: WindramI 
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119681 : iMilev et all Il975| ). High-resolution observations 
of the head structure of NGC 1265 reveal two radio 
jets emerging from the galaxy which are at an angle 
of 90° to the tail where the tw o jets are flrst paralle l 
to each other a nd th en merge (jO'Dea fc OwenI |1986[) . 
iGisler fc MilevI ()1979() found an extension of this ra- 
dio tail to the north-east. Surprisingly, a t lower sur- 
face brightness, iSijbring fc de Bruvnl (|1998[ ) flnd an ad- 
ditional large scale structure that arches toward the east 
around the steep spectrum tail of NGC 1265 and forms 
a closed ring with a very steep, but constant spectral 
index across. Using rotatio n measure (RM) synthesis, 
Ide Bruvn fc BrentiensI ([20051 ) flnd a very weak diffuse po- 
larized structure with a Faraday depth of approximately 
50 radm~^ at the angular position coincident with the 
steep spectrum tail as well as the ring structure of NGC 
1265. Follow-up mosaic observations of the area around 
the Perseus cluster show t hat most if not all of this emis- 
si on is of Galactic origi n (li3rentien£ll201lD . 

iSiibring fc de Bruvnl (| 19981 ) discuss a total of four 
models to explain this puzzling observation. They im- 
mediately discard two models — one model of chance su- 
perposition of several independent head-tail galaxies due 
to the lack of strong radio sources in this field (Model 
1) as well as another model that hypothesizes reaccel- 
eration of mildly relativistic electrons in the turbulent 
wake of a galaxy due to contrived projection probabilities 
and implausible energetics (Model 2). The previously 
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favore d pair of models discussed in iSiibring fc de Bruvnl 
(|1998D — apparently simple — do in fact appear to have 
considerable amount of fine-tuning and have severe phys- 
ical problems associated with them as we will point out 
here. Model 3: in this model, they postulate a heli- 
cal cluster wind that has to be aligned with the line-of- 
sight (LOS) to produce the observed ring structure in 
projection: such a configuration is highly unlikely, and 
this model has problems in explaining the constancy of 
the spectrum and the surface brightness along the radio 
ring. To balance the synchrotron and inverse Compton 
cooling of an electron population, this model would need 
to postulate a fine-tuned acceleration process that how- 
ever must not fan out the well-confined radio emission 
along the arc. Model 4: in this model, the "radio tail" 
would outline the ballistic orbit of NGC 1265. To sus- 
tain such a helical orbit of NGC 1265 over 360°, this 
model would require an undetected dark object of mass 
M > Mngc 1265 =i 3 X W^Mq orbiting the galaxy.^ 
The change of the direction in the brighter part of the 
tail remains unexplainable and there is also the prob- 
lem in explaining the constancy of the spectrum and the 
surface brightness along the radio ring. Given these dif- 
ficulties, it is attractive to consider alternative possibil- 
ities that may explore the interaction between a radio 
galaxy with the outskirts of the Perseus ICM. This has 
been f oreshadowed in a remark bv lde Bruvn fc BrentiensI 
(|2005D who speculate whether the large scale polarized 
structure that arches around the steep spectrum tail of 
NGC 1265 is indeed the remains of an earlier phase of 
feedback. Our work will demonstrate that this picture 
is not only the simplest consistent explanation for the 
radio morphology and spectrum, but we also use it to 
indirectly infer the presence of a cluster shock wave and 
measure its properties. 

Previously, the morphology of a giant radio galaxy 
has already been used to indirectly detect a large scale 
shock at an intersecting filament of galaxies (lEnfilin et al.l 
I2001D by using the radio g alaxy as a giant cluster 
weather station ([Burnsl [l998[ ). Gravitationally driven, 
supersonic flows of intergalactic gas follow these fila- 
ments toward clusters o f galaxies that r epresent the 
knots of the cosmic web (jBond et al.lll996h . The flows 
will inevitabl y collide and form large-scale shock waves 
(lOuilis et al.iri 998: Mi niati et al1[2000t iRvu et al.ll200l 
iPfrommer et al.i ,2006. ) . Of great interest is the sub- 
class of accretion shocks that are thought to heat the 
baryons of the warm-hot intergalactic medium (IGM) 
when they are accreted onto a galaxy cluster. Forma- 
tion shocks have also been proposed as poss ible genera- 
tion sites of inte rgalact i c mag netic fields (|Kulsrud et al.l 
Il997t IRvu et al.lll998al 120081) . Prior to this work, only 
discrete mer ger shock waves have b een detected in the 
X-rays fe.g.. iMarkevitch et al.l[2002| ) and there was no 
characterization possible of the detailed flow properties 
in the post-shock regime as to test whether shear flows — 
the necessary condition for generating magnetic fields — 

^ Since the velocity dispersion of NGC 1265 is unknown, we use 
the Faber-Jackson rela tion in combination with the My,^-a relation 
(|Tremaine et al.|[2003 ) to obtain the stellar mass of NGC 1265 of 
M* ~ 2.5 X 10"'^"'^ Af0. We estimate the total mass by using the 
universal baryon fraction and a stellar-to-baryon mass yield of 0.5 
so that the resulting halos mass should be considered a lower limit 
to the true one. 



are present. 

The structure of the paper is as follows. In Section [21 
we present the basic picture of our model in a nutshell 
while we derive the detailed three-dimensional (3D) ge- 
ometry of NGC 1265 within our model in Section [3] In 
Section 21 we work out the properties of the accretion 
shock onto the Perseus cluster including those of the 
post-shock fiow and present the implications for the in- 
falling warm-hot IGM. In Section |S1 we carefully exam- 
ine the underlying physics as well as the hydrodynamic 
stability of our model and discuss our findings in Sec- 
tion |6l Throughout this work, we use a Hubble constant 
of Hq — 70km s^^Mpc^^. For the currently favored 
ACDM cosmology with the present day density of to- 
tal matter, ilm = 0.28, and the cosmological constant, 
J7a — 0.72, we obtain an angular diameter distance to 
Perseus {z = 0.0179) of Dang = 75Mpc; at this distance, 
1' corresponds to 21.8 kpc. X-ray data estimates a virial 
radius and mass for Per seus"^ of i?2on = 1 -9 Mpc and 
= 7.7 X 10"Afo fReiprich fc Boliringed[200a) . 

2. THE IDEA IN A NUTSHELL 

Before we present the idea of our model, we point out 
the main morphological and spectral properties of the 
giant radio galaxy NGC 1265 that every model would 
have to explain: (1) The synchrotron surface brightness, 
S'tx, and the spectral index, a, between 49 and 92 cm 
along the tail of NGC 1265 (starting at the galaxy's head) 
sh ow a characteristic behavio ur (as shown in Figure 2 
of ISiibring fc de Bruvnl H 9981 ). In the first part of the 
tail, both quantities decline moderately in a way that is 
consistent with synchrotron and inverse Compton cooling 
of a relativistic electron population that got accelerated 
at the base or the inner regions of the jet. (2) At the point 
of the tail where the twist changes in projection from 
left- to right-handed, both quantities experience a sudden 
drop — while 5^ changes by a factor of 10, a declines from 
— 1.1 to —2.1. (3) Finally, along the remaining curved 
arc, Sy and a stay approximately constant on a total 
arc length of Zarc = 27ri?^ vT+A? > 700 kpc, where R ~ 
150 kpc is the radius of the arc, ^ ~ 3/4 the projected 
arc length in units of 27r radians, and k = h/ {2ttR) > 0, 
where h is the height of the 3D helix. This property is 
in particular puzzling, as there is no visible synchrotron 
cooling or fanning out of the dilute part of the tail visible. 

These findings taken together suggest the presence of 
two separate populations of relativistic electrons giving 
rise to the bright and the dim part of the tail, respec- 
tively, where the latter must have experienced a coherent 
energetization event^ over a length scale of 2R ~ 300 kpc 
and on a timescale that is shorter than the cooling time of 
the radio emitting electrons of Tgync, ic ^ 2.9 x 10® yr. The 
presence of one radio tail that connects the two electron 
populations in projection points to a causally connected 
origin of the synchrotron radiating structure. The most 

We define the virial mass A/^ and the virial radius /J^ as 
the mass and radius of a sphere enclosing a mean density that is 
A = 200 times the critical density of the Universe. 

^ We discard the alternative that there is a continuous accelera- 
tion process that exactly balances synchrotron and inverse Comp- 
ton cooling due to the large degree of fine-tuning needed and be- 
cause it is difficult to avoid a fanning out of the synchrotron emis- 
sion region along the tail if the acceleration process is, e.g., of 
turbulent nature. 
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Fig. 1. — Schematic snapshots of the time evolution of our model of NGC 1265 (cut plane of the top view with the two axes representing 
the LOS and a suitably chosen angular direction, not to scale). Frame shows the initial configuration of NGC 1265 together with two 
detached plasma bubbles from previous outbursts as they are moving upward toward the shock surface. In frame 1, the first bubble on 
the right passes through the shock and transforms into a torus (vortex ring) which re-energizes an aged electron population by means of 
adiabatic compression so that it emits observable radio synchrotron radiation. As shown in frame 2, its orbit is subject to the coherent 
shock deflection upon oblique shock crossing. In this frame, NGC 1265 also crosses the shock but remains on its ballistic orbit. However 
in our model, the shock passage triggers a new outflow that is now seen as head-tail radio emission. The jets are being bent by the ram 
pressure wind which causes them to merge into a single radio-emitting tail at low radio frequencies (shown in the schematic inlay in frame 
3). Also shown in this frame is the transformation of the second plasma bubble on the left into a torus. Frame 4 shows today's configuration: 
the head-tail jet of NGC 1265 is shaped by post-shock flows owing to the oblique shock crossing, the second torus on the left is stretched 
by post-shock shear flows (Figure [Sjl and emits diffuse soft radio emission, and the electron population in the first torus on the right (not 
drawn) has sufficiently cooled such that its radio emission is too weak to be detected by current radio telescopes. 



natural explanation that combines these observational 
requirements are the reminders of two distinct epochs of 
active galactic nucleus outbursts where the most recent 
one is still visible as a head-tail radio jet and the older 
one experienced a recent coherent energetization event. 
In particular, we propose that such an energizing event 
could be provided by the passage of a detached radio 
plasma bubble from a previous outburst through a shock 
wave. This passage transforms the plasma bubble into 
a torus (vortex ring) and adiabatically compresses and 
energizes the aged electron population to emit low sur- 
face brightness and steep-spectrum radio emission (as we 
detail below). If the shock crossing is oblique, then the 
radio torus and the tail end of the recent outburst ex- 
perience the same degree of shock deflection since both 
structures can be regarded as passive tracers of the post- 
shock velocity field and hence provide a natural expla- 
nation for the apparent connectivity of both structures. 
With this idea in mind, we sketch a schematic of the time 
evolution in our model in Figure [TJ 

The reason for the transformation of a plasma bub- 
ble into a toroidal vortex ring can be easiest seen in the 
rest frame of the shock, where the ram pressure of the 
pre-shock gas balances the thermal pressure in the post- 
shock regime. The bubble is filled with hot (relativistic) 
and more dilute plasma compared to the surrounding 
ICM. Once the dilute radio plasma of the bubble comes 
into contact with the shock surface, the ram pressure 
is reduced at this point of contact due to the smaller 
density inside the bubble (EnBlin fc B riiggen 2002) . The 
shock and the post-shock gas expand into the bubble 
and propagate with a faster velocity compared to the 
incident shock in the ICM. Owing to symmetry, the am- 
bient gas penetrates the line through the center of the 
bubble first and has a smaller velocity for larger im- 
pact parameters. This difference in propagation veloc- 
ities implies a shear flow and eventually causes a vortex 
flow around the newly formed torus which stabilizes it 
as it moves now with the post-shock velocity field. To 



study the timescale on which this transformation hap- 
pens, we solve the one-dimensional (ID) Riemann prob- 
lem of a shock passage through a bubble exactly in Ap- 
pendix [BJ These calculations are complemented by a 
suite of two-dimensional (2D) axisymmetric simulations 
using a code that employs an upwind, total variation 
diminishing sch eme to solve the h ydrodynamical equa- 
tions of motion (|Rvu et al.lll998bD . To map out param- 
eter space, we varied the Mach number and initial bub- 
ble/ICM density contrast (see Figure [2] for one realiza- 
tion). As a result, we find that an initially spheroidal 
bubble will then evolve into a torus on a timescale Tform 
that is determined by the crossing time of the origi- 
nal bubble-ICM contact discontinuity (CD) through the 
bubb le. F or typical numbers, Tform — 1.4 x 10* yr (Equa- 
tion (IBGI) ). As we will see in Section[5l this is much faster 
than any transverse shear on the scale of the bubble or 
beyond can act to distort it. 

The morphology of the radio torus is consistent with 
MHD simulations of t his effect that assurn e strong shocks 
(cf. Figs. 9 and 10 in lEnfilin fc Bruggen| [2002). The rel- 
ativistic plasma within the bubble experiences an adi- 
abatic compression by a factor^ C that increases the 
Lorentz factor of the relativistic electrons as F c>c C^/^ 
and the rms magnetic field as i? esc C^/'^. Hence the 
radio cutoff of a cooled electron population increases as 
^'max oc BT^ cx C*/^ and illuminates a previously un- 
observable radio plasma (Enfilin &__Gopal-Krishna 2001). 
Synchrotron and inverse Compton aging develops a steep 
spectrum with spectral index a ~ 2 and a low surface 
brightness Si, cx v^"^. 

The galaxy itself remains on a ballistic orbit that is 

^ In this work we will encounter three different compression fac- 
tors/density contrasts that should not be confused: the compres- 
sion factor C = Vbubblc/Horus by which the volume of the bubble 
changes upon shock passage, the compression factor C'b = P2/P1 
between the post- and pre-shock ICM density at the shock front, 
and the density contrast 5 = ntorus/'^icm between the density of 
the radio plasma and the ICM. 
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Fig. 2. — Left two panels: images from a 2D axisymmetric simulation showing the evolution of a spherical bubble overrun by a Mach 5 
shock that transforms into a torus (each half-panel shows a selected snap shot). Base 10 log of gas density is shown. The initial bubble/ICM 
density contrast <5 = 10~^, the bubble diameter is unit length, the pre-shock ICM density is unity, as is the pre-shock sound speed. Times 
correspond to initial shock impact on the bubble (t = 0), just before the internal bubble shock has crossed the bubble (not visible with 
this color scale, t = 0.025), the ICM shock has crossed a bubble radius (externally, t = 0.1) and just after the bubble contact discontinuity 
has crossed the bubble {t = 0.175). Right panel: the velocity field in the rest frame of the ICM post-shock fiow overlaid by gas density 
contours. Note that the stabilizing vortex fiow around the torus has already established at the time when the torus has formed. 



unaffected by the shock passage that fifcely triggers a 
new outflow that is now seen as a head-tail radio emis- 
sion with the jets being bent by the ram pressure wind. 
The associated compression wave propagating through 
the interstellar medium could have triggered a bar-like 
instability in the inner accretion disk that enabled effi- 
cient angular momentum transport outward and accre- 
tion onto the super-massive black hole (SMBH) which 
was responsible for launching the jet. 

This model explains the observed spectral steepening 
of the head-tail radio galaxy of Aa ~ 0.5 due to radia- 
tive cooling along the bright part of the tail. In par- 
ticular, our model naturally accounts for the observed 
sudden steepening of the spectral index and dropping 
synchrotron brightness at the transition (in projection) 
of the radio tail of the current outburst to the shock- 
illuminated radio torus and match the observed con- 
stancy of spectrum and surface brightness across the 
torus (see Figure [3] and Siibring & de Bruyn 1998). 

The shock compression acts mostly perpendicular to 
the filaments during the formation of the torus. Con- 
sequently, the field component parallel to the bubble 
surface is preferentially amplified an d dominates even- 
tually the magnetic energy density (|EnBlin fc Briiggenl 
12002!) ■ (Provided the initial magnetic field component 
parallel to the shock was not too weak to be ampli- 
fied above the perpendicular components.) Hence, our 
model predicts radial synchrotron polarization vectors 
relative to the center of the radio torus with a polar- 
ization degree of ~ 5% for our derived geometry of a 
viewing angle of x ~ 23° (see Section [3]). This is consis- 
tent with upper limits on the polarization of 10%'-25% 
(jSiibring fc de Bruvnlll998l ). Future low- frequency radio 
observations with, e.g., LOFAR should be able to detect 
this signal. 

3. OBSERVABLES AND GEOMETRY 

In this slightly technical section, we derive the spe- 
cific geometry that we will use later on to demonstrate 
the physical plausibility of our proposed model — if the 
reader is only interested in the physical ideas and re- 



sults, this section might be easily skipped without loss of 
the underlying logic of this work. 

In general, the location and orientation of the shock 
surface with respect to the galaxy's infall velocity vec- 
tor and the projection geometry of this system are de- 
generate. Inspired by cosmological simulations that 
quantify structure formation shocks (iRvu et al.l 120031 : 
iPfrommer et all 120081 : iHoeft et al.l |2008D, we take the 
plausible assumptions that (1) the shock surface is 
aligned with the gravitational equipotential surface of the 
Perseus cluster that we determine by taking concentric 
spheres around NGC 1275 and (2) the second outburst 
giving rise to the currently observed head-tail structure 
has started shortly after shock crossing so that we can 
identify both events. 

3.1. Observables and Simple Derived Quantities 
3.1.1. The Perseus accretion stiock 

NGC 1265 has a large radial infall velocity of Vr = 
2170 km s^^ that has been obtained after subtracting 
the mean heliocentric velocity of the Perseus cluster of 
Vr,Per — 5366km s"^. The Faraday RM across th e cen- 
tral few kpc of the two jets (|0'Dea fc Owenlll986D show 
a scatter of 20 radm~ around a mean value of around 
25 radm"^. The scatter is most probably due to the 
magnetized plasma in the cocoon surrou nding the jet or 
the i nterstellar medium of NGC 1265 (jO'Dea fc OwenI 
Il987f ). The low value of RM fluctuations strongly ar- 
gues for an infalling geometry with NGC 1265 sitting 
on the near side of Perseus and little Faraday rotating 
material in between us and the source. This argument 
strongly favors the picture that the shock that trans- 
formed the radio bubble is the accretion shock rather 
than a merger shock. The strongest argument in favor 
of this picture would be the identification of the filament 
along which NGC 1265 was accreted in redshift space. 
Since NGC 1265 is only 600 kpc offset from the center of 
the Perseus cluster (in projection), the finger-of-God ef- 
fect makes it very challenging to isolate a galaxy fllament. 
There is however indirect evidence for such a cosmic fila- 
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TABLE 1 
Definitions 



er, = 1 LOS vector pointing to observer 

Us, nj. = 1 Shock normal, pointing outward, away from convexly curved surface 

V, = 3D proper motion of NGC 1265 relative to Perseus, pointing inward 

X = arccos( e,. ■ Us) "Shock orientation," angle between LOS and shock normal 

9 = arccos(— er ■ ^) "Inclination of the galaxy's orbit," angle between negative LOS and galaxy velocity 

(j> = arccos(— ris ■ — ) "Shock obliquity," angle between negative shock normal and galaxy velocity 

Subscripts t, r Transverse (on the plane of the sky), radial (along LOS) 

Subscripts ±, || Perpendicular, parallel to shock normal 

Subscripts 1, 2 Pre-shock, post-shock regime; refers to quantities of the gas 



ment from Perseus X-ray data. Dark Ha filaments near 
the LOS to the core do not show any foreground cluster 
X-ray emission and have redshifted velocities of around 
3000 km s~^ relative to the Perseus cluster, so they ap- 
pear to be falling in, too (A. Fabian, private communi- 
cation) . 

3.1.2. Electron cooling timescales 

The radio synchrotron radiating electrons of Lorentz 
factor r emit at a frequency 



^syn 



3eSr^ 

2'KmeC ' 



(1) 



where e is the elementary charge, rrie is the electron mass, 
and c is the light speed. Synchrotron and inverse Comp- 
ton aging of relativistic electrons occurs on a timescale 



'syn, ic — / „2 



cmb 



s2)r' 



(2) 



where -Bcmb — 3.2/^G (1 -I- z)^ and ctt is the Thompson 
cross section. Combining both equations by eliminating 
the Lorentz factor F yields the cooling time of electrons 
that emit at frequency Vsym 



2 

cmb 



<2.9xl0*^yr. (3) 



The highest frequency I'syn = 600 MHz at which the torus 
can be observed gives the shortest cooling time for any 
given magnetic field value. Interestingly, Tgyn, ic is then 
bound from above and attains its maximum cooling time 
'At B = Bcmb/\/3 ~ 1.8 txG (1 + zf. 

3.1.3. Ltmits on the galaxy's inclination angle 

We denote the projected length on the plane of the 
sky that the galaxy has traveled since shock crossing 
by Lf gal ■ Looking at the bending of the jet of NGC 
1265 (jo Pe a fc O wenlll986D and the o verall morphology 
of the tail (Siibrin g fc de Bruvnlll998l) , we conclude that 
the past orbit of NGC 1265 in the plane of the sky is 
directed mostly southward. We obtain a lower limit on 
Lt,ga\ ^ 220 kpc by measuring the vertical distance of the 
"h ead" to the end of the bri ght radio tail (Figure 1(a) 
in iSiibring fc de Bruvnl [19981) . Assuming that the time 
of shock passage of the galaxy and the radio bubble co- 
incide, we can obtain a lower limit on the inclination of 
the galaxy's with the LOS, 



9 > arctan 



min(Lt,gai) 



Vj. max(Tsyn^ ic) 



19°, 



(4) 



where Vr = 2170 km s^ . We will see that detailed geo- 
metric considerations of the overall morphology of NGC 
1265 show that there is a time lag between the shock 
passage of the galaxy and the bubble of r ~ 6 x 10^ yr 
that one has to add to Tgyn^ ic such that the limit on the 
inclination weakens slightly and becomes 9 > 16°. 

3.2. Geometrical Constraints 

It turns out that the morphological richness of NGC 
1265 in combination with physical arguments then con- 
strain the geometry and shock properties of this system 
surprisingly well. To this end, we derive a general expres- 
sion for the gas velocity after shock passage by taking 
into account the shock obliquity, 0, the shock orienta- 
tion, X, and the inclination of the galaxy's orbit, 9 (see 
Table[l]for definitions). Using mass and momentum con- 
servation at an oblique shock with a compression factor 
Cs = P2/P1, the velocities parallel and perpendicular to 
the shock read as 



V + V COS (p Tls 

V 



V2, 



V2 — V+V COS ( 



Cs 



(5) 
(6) 

(7) 



In order to connect to observables, we derive the velocity 
components parallel and perpendicular to the LOS, 



V2,r ~ Sr ■ V2 = Vr { I 



V2,t 



cos X COS Cs 



1 



{V2 - V2,rery 



COS f 
0,2 



COS^ 



o / COS X COS (f) Cs ~ 1 



COS 9 



Cs 



1 — COS^ 
21 1/2 



(8) 



C2 



(9) 



We will determine the shock compression factor Cs in 
Section|4]based solely on the volume change of the bubble 
while assuming that the bubble had negligible ellipticity 
prior to shock crossing. 

We employ the following strategy to determine the 
geometry. (1) We choose the accretion shock radius 
Rs ^ -R200 — 1.9 Mpc since the hot gas within the clus- 
ter is unlikely to support a shock with Mach number 
~ 4. (2) A priori, the position of the shock pas- 
sage of NGC 1265 and its advected radio bubble is un- 
known. We choose the point where the galaxy passed 
the shock as the origin of the coordinate system O, use 
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the current galaxy's position as a starting value, and in- 
terpolate the past orbits of the galaxy and the tail in 
the plane of the sky. An iteration of the position of O 
(that involves a complete cycle through the points listed 
here) quickly converges on a consistent model. Subject 
to our assumption of the shock surface, this yields the 
shock normal, rig, the shock orientation, x, and projected 
orbit that the galaxy has traveled since shock crossing, 
it, gal- (3) Choosing the inclination of the galaxy's or- 
bit, 9, determines the 3D velocity vector v and the time 
since shock crossing, = Lt^gai/(urtan0). (4) Momen- 
tum conservation at the shock implies for the post-shock 
gas that its velocity vector lies in the plane that con- 
tains ris and v and is defined by r ■ (n^ x — 0. 
The intersection of this plane and the plane of the sky 
(given by z = for our choice of coordinates) yields 
the past transverse orbit of the post-shock gas, V2,t- (5) 
Solving Equation Q for the unknown shock obliquity 
(f> enables us to derive the current radial position of 
the galaxy (relative to the cluster center) according to 

''gal — Rs — VTs cos (f> with V = VrVl + tau^ 9. (6) We 
compare the galaxy's velocity, v, to the escape velocity, 
Wosc(?'gai) = ^/2GM{<r^^l)Jr^l. For the large radial ve- 
locity of Vr = 2170 km s^^, it is not trivial to meet the 
criterion, v < t'esc, stating that the galaxy is gravita- 
tionally bound; hence we prefer small values of Rg and 
9 < 32° (while still evading the electron cooling bound 
on 9 as derived in Equation (jl])). The head-tail morphol- 
ogy of the jet also argues for a large ram pressure that 
it can only experience inside R200 and hence for small 
values of Rg. We finally check whether the model vio- 
lates any constraints on the shock obliquity (j) as will be 
derived in Section 15.21 If this does not yield a consistent 
model we vary the shock radius Rs , the inclination of the 
galaxy's orbit 9, and the position of the galaxy's shock 
crossing, O, and start the next iteration until we arrive 
at a self-consistent and physically plausible model. 

Our final model that fits best these constraints is 
shown in Figure |3] and has the parameters Cg = 3.4, 
Rg = i?2oo = 1.9 Mpc, 9 = 32°, = 9°, X = 23°. 
We find the time since shock crossing of NGC 1265 to 
be 1.8 X 10* yr, the Cartesian vector of the shock nor- 
mal n = (0.2,0.34,0.92), and the velocity vector for 
NGC 1265 oi V = (O,sin6',cos6') (note that we de- 
fine the Cartesian coordinate system in Figure ^ . The 
galaxy's velocity has a radial and transverse compo- 
nent of Vr = 2170km s~^ and Vt = 1360km s~^, re- 
spectively, which yields a total velocity of the galaxy, 
V = 2550 km s~^. This is only slightly larger than the 
escape velocity, Wcsc(?'gai) — 2350km s~^, where we used 
rgai = 1.45 Mpc, M(< rgai) ~ XturbM200 - 9 X 10"Mq, 
neglect a logarithmic correction factor of the mass and 
assume a turbu lent pressure support of Xturh = 0. 2 
(|Rvu et al.ll200l iLau et all 120091 : iBattaglia et al.llMoh . 
Note that if one relaxes the escape velocity constraint 
and assumes that tidal processes are able to dissipate 
more energy to bind the galaxy during first passage, we 
find consistent solutions at larger shock radii that nev- 
ertheless lie close-by in the parameter space {9,(j),x,0). 
The choice of the virial radius as the site of the accre- 
tion shock might seem to be too small compared to the 
typical locations inferred from cosmological simulations 
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Fig. 3. — Plane of the sky view (top) and view from above (bot- 
tom) of the time evolution of NGC 1265 through the accretion 
shock of Perseus (radio su rface brightness at 600 MHz taken from 
ISiibri ng fc de Bruvnlll99SI ) . Note that the angles of the cut-planes 
through the shock surfaces (at the point of shock passage of NGC 
1265, labeled with 2) and velocity vectors are derived from our 
best-fitting geometry, are drawn to scale (given at the top right), 
and the numbers correspond to the time sequence defined in Fig- 
ure [1] The galaxy and the detached radio bubble are advected 
in a filament (initial configuration labeled with 0, bubble on the 
left shown in projection as orange half-circle, bubble on the right 
in Figure [1] omitted for clarity). The galaxy stays on its ballistic 
orbit unaffected by the shock passage (2). In the top panel, the 
shock surface is almost face-on with a small inclination of x = 23° 
to the LOS whereas the bubble is offset from the cut plane and sit- 
uated closer to the observer. This inclination of the jet axis from 
the plane of the sky implies a later shock crossing time of the bub- 
ble compared to the galaxy. After T2_>3 ~ 6 x 10^ yr, the bubble 
passes through the shock and transforms into a torus (3). The ra- 
dio torus and the tail of NGC 1265 experience an almost coherent 
deflection due to the oblique shock passage for the remaining time 
of T3_j.4 ~ 1.2 X 10* yr until today (4). The red circle shows the 
bubble circumference when the shock surface comes into contact 
with the bubble (3, shown as light red point); due to the inclina- 
tion of the shock with respect to the LOS, this point would not be 
visible in projection as it appears in both cases at the back side of 
the bubble. The fact that projected shock surface coincides with 
the rim of the torus is accidental and depends on the choice of the 
radial position of the shock Rs which we choose as Rs ^ R'200- 
A shell of 250 kpc in the region interior to the accretion shock 
appears to experience an additional shear flow caused by vorticity 
generation at the curved shock surface. 
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(iMiniati et aLl 120001 : iKang et aLl 120071 : iPfrommer et all 
12008'). However, the location of the accretion shock (in 
particular along a filament) is not stationary but dynam- 
ically determined. Depending on the ram pressure of the 
accreting material and the post-shock pressure, its posi- 
tion will re-adjust dynamically to account for the conser- 
vation laws. We have shown that the velocity of NGC 
1265 is quite large which implies a large ram pressure 
and hence should yield a shock position that lies closer 
toward the cluster center compared to the shock position 
where dilute IGM accretes from voids with a considerable 
smaller density and ram pressure. We note that there 
have been attempts to constrain the 3D velocity of NGC 
1265 based on the observed sharpness of its X-ray surface 
brightness map and the optical v elocity dispersion of the 
Perseus cluster (|Sun etall [20051) . To arrive at a lower 
bound on the orbit inclination, 6'x-ray > 45°, these au- 
thors use an axially symmetric density profile that how- 
ever does not take into account magnetic draping. This 
naturally would account for the sharpened interface be- 
tween the ICM and the interstellar medium and might 
break the assumed symmetry depend ing on the morphol- 
ogy of the ambient mag netic field (|Dursi fc: Pfromrnerl 
mm IPfrommer fe Dursill2ninn . 

The post-shock gas velocity — as traced by the radio 
emitting tail — is reduced at the shock by transferring a 
fraction of the kinetic to internal energy. For our choice 
of parameters, we find V2 = w^/l^^^os^^TZTj^^lJ/Cf ~ 
850 km s~^ that divides into the radial and transverse 
LOS component of the tail as V2,r — 520km s~^ and 
U2,t — 670 km s~^, respectively. The velocity compo- 
nents parallel and perpendicular to the shock normal 
are Vi^\\ = vcoscj) ~ 2530km s^"'^, W2,|| — vcos(f)/Cs — 
740km s~^, and v±_ — vi^j_ — W2,_l — vsinc/) ~ 
400 km s~^. Once in the post-shock regime, the radio 
emitting tail and the torus experience the same large- 
scale velocity field. To explain the observed direction of 
the bending of the tail toward west, the shock normal 
needs to have a component pointing westward (for the 
observed direction of the galaxy's velocity southward). 

What about the orientation of the line connecting the 
bubble center and NGC 1265 (shown as dashed lines in 
Figure [S])? Assuming that the spin orientation of the 
SMBH as possibly traced by the radio jets did not change 
during shock passage, the previous jet blowing the radio 
bubbles predominantly had an east-west (EW) compo- 
nent in the plane of the sky. If the jet angle had no com- 
ponent along the LOS, the bubble had crossed the shock 
earlier than NGC 1265 due to the assumed convex shock 
curvature. Hence the bubble would need to have been at 
a larger distance from NGC 1265 than in our model and 
would have needed a longer time r ~ 3 x 10^ yr > Tsyu, ic 
to reach the current position (extrapolating the position 
back to the shock with the orientation of the deflection 
vector given by our model). It follows that the shock- 
enhanced radio emission would have faded away by now 
and the bubble would have had to spend at least a factor 
of 2.5 longer in the post-shock region which would not be 
consistent with the assumed EW alignment of the previ- 
ous jet axis. Hence our model postulates that the eastern 
bubble was on the front side relative to the observer and 
had a smaller LOS than NGC 1265 (see Figure[3l). This 



caused a "two-stage process": the bubble continued for 
T2^3 ~ 6 X 10^ yr on its orbit in the filament after NGC 
1265 passed the shock until it came into contact with the 
shock surface. This configuration has the advantage of 
the faster pre-shock velocity and associated larger length 
scale in the plane of the sky that the bubble traversed 
and enables us to fulfill the geometric requirement to 
cross the shock surface while ending up at current posi- 
tion. After shock passage and its transformation into a 
torus, it experiences the an almost similar deflection as 
the tail of NGC 1265 due to the oblique shock passage 
for the remaining time of T3-^4 ~ 1.2 x 10® yr until today. 
This also solves the problem of the non-observation of 
the "second" bubble that the previous western jet should 
have blown: in our proposed geometry, this hypothetical 
bubble was situated westward of NGC 1265 at the far 
side of the cut plane in Figure |3l It should have passed 
the shock r ^ 2 x 10* yr earlier than the eastern bub- 
ble and faded away by now (Equation ([3|)). Note that 
the position and morphology of the accretion shock is a 
function of time: as the radio bubble comes into contact 
with the shock surface, the shock expands into the light 
radio plasma toward the observer and attains additional 
curvature. For simplicity we suppress this effect in our 
sketch in Figure [3] and place the torus in the lower panel 
at a distance from the (initial and unmodified) shock sur- 
face which should be equal to that of the modified shock 
surface. 

4. PROPERTIES OF THE PERSEUS ACCRETION 
SHOCK 

Using the 49 cm image of NGC 1265 (Figure 1(a) in 
ISiibring fc de Bruvnl [TQQI . we measure the major and 
minor radius of the torus to i? ~ 7' and rmin — l.'l — 1.'4. 
We only consider orientations that lie in the northern and 
southern sectors as the elongated appearance toward the 
EW might be caused by shear flows as we will argue 
later on in this section. At the distance of Perseus, this 
corresponds to physical length scales of i? ~ 150 kpc and 
r,„iii ~ (24 . . . 30) kpc, respectively. Assuming that the 
bubble was spherical before shock crossing and that the 
major radius did not change (Enfilin & Briiggen 2002), 
we estimate a compression factor through the associated 
volume change of 

V^bubbie ^ l^rR' 2 ^ 

^miii / 

The radio plasma is adiabatically compressed across the 
shock passage according to P2 /Pi — C^'"^ , where 7roi — 
4/3 for an ultra- relativistic equation of state. However, 
the adiabatic index can in principle approach 7^01 ^5/3 
if the bubble pressure is dominated by a hot but non- 
relativistic ionic component. Varying 7101 — 1.33 — 1.5 
and C ~ 6 - 10, we obtain P2/P1 ~ 21.5 ± 10.5. As- 
suming that the radio bubble is in pressure equilibrium 
with its surroundings before and after shock crossing, 
this pressure jump corresponds to the jump at the shock. 
Applying standard Rankine Hugoniot jump conditions 
(jLandau fc Lifshitzlll9591) of an ideal fluid of adiabatic in- 
dex 7 = 5/3 yields a shock Mach number of ~ 4.2+°;2, 
a density jump of 3.4lg'4, and a temperature jump of 
6.3^^2 J- We assumed a flat prior on the uncertainties of 
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C and 7ici which yields symmetric error bars around the 
pressure jump. We quote the jumps in the other ther- 
modynamic quantities corresponding to our mean PijPx- 
These quantities inherit asymmetric error bars due to 
the non-hnear dependence of these jumps on Mach num- 
ber. We caution tha t we performed these estima tes from 
the 600 MHz map (jSiibring fc de Bruvn|[T99l ): future 
higher resolution and more sensitive radio observations 
are needed to assess the uncertainties associated with fi- 
nite beam width and to obtain a more reliable morphol- 
ogy of the low surface brightness regions of the torus. 

Extrapolating X- ray profiles of Perseus 
()Churazov et al.l l2003| ) to i?2oo = 1-9 Mpc and us- 
ing our mean values for the jumps in thermodynamic 
quantities, we can derive pre-shock values for the 
gas temperature and density that reflect upper lim- 
its on the gas prop erties in the infalling warm-hot 
IGM as predicted by iBattagha et"all (|2009r ). We find 
fcTi < 0.4 keV, ni = 1.93 < 5 x IQ-^cm-^, and 
P\ S 3-6 X lO"^** ergcm"'^. For the cluster temperature 
profile, we combine the central profile from X-ray 
data with the temperature profiles toward the cluster 
perip hery according to cosmological cluster simulation s 
(e.g., iPfrommer et al.1 120071: iPinzke fc Pfrommeii [2010l ) . 
yielding 



fcTn 



kTo 



1 



0.2 R 



200 



-0.3 



l + (r/rc)- 

' (11) 

where kTo — 3keV, fcTmax = 7keV, and rc = 94kpc. 
This represents the correct temperature profile over the 
entire range of the cluster when comparing to mosaiced 
observations by Suzaku (S. Allen, private communica- 
tion). 

We now present evidence for post-shock shear flows 
with an argument that is based on the projected orien- 
tation of the ellipsoidally shaped radio torus whose main 
axis is almost aligned with the EW direction. We con- 
sider two cases. (1) If the ellipticity was due to projection 
of a ring-like torus, the shock normal Us could not have 
components in the EW direction. Since the transverse 
velocity of the galaxy is pointing southward, momentum 
conservation at the (oblique) shock would imply a deflec- 
tion of the post-shock gas in the plane that contains 
and V — without components in the EW direction. Hence 
an additional shear flow would be needed to explain the 
westward bending of the galaxy's tail even though the 
smooth and coherent bending of the tail might be diffi- 
cult to reconcile with the turbulent nature of vorticity- 
induced shear flows. (2) If the bending of the tail was due 
to oblique shock defiection, Ug would have to have a com- 
ponent pointing westward. Projecting an intrinsically 
ring-like torus — aligned with the shock surface — would 
yield an apparent ellipsoidal torus with the main axis at 
some angle with the EW direction on the plane of the sky. 
This argues for a shear flow that realigns the orientation 
of the ellipsoid with the observed EW direction. In fact, 
assuming that the shock surface is aligned with the grav- 
itational equipotential surface of the Perseus cluster, one 
can show that the implied curvature of the shock surface 
causes a vorticity in the post-shock regime that shears 



the post-shock gas westward^. We cannot completely 
exclude pre-existing ellipticity in the radio bubble but 
flnd it very unlikely to be the dominant source of the ob- 
served ellipticity as the same directional deflection of the 
tail structure as well as the radio torus would then be a 
pure coincidence rather than a prediction of the model. 

The observed flattening of the projected radio torus 
amounts to / = 1 — 6/a ~ 0.3, where a and b are the 
major and minor axes of the ellipse, respectively. With 
the assumption that the shock surface is aligned with 
the gravitational equipotential surface and using the re- 
sult from Section [31 the flattening due to projection, 
f — 1 — coacf) ~ 0.066, falls short of the observed one 
calling again for shear flows to account for the differ- 
ence. The amount of vorticity uj injected at a curved 
shock of curvature radius Rcmv should be proportional 
to the perpendicular v elocity v± over the bubble radius 
^'bubble (ILighthilllllQSTi) . 



Rc 



= /- 



^'bubble 



(12) 



where the constant of proportionality / depends on the 
flattening / of the torus. Using / ~ / ~ 0.3, Cs — 3.4, 
and Tbubbio = 150 kpc, we obtain a rough estimate for 
the curvature radius of the shock of i?curv — 850 kpc 
which is also adopted in Figure [31 This is similar to the 
dimensions of a filament connecting to a galaxy cluster 
and provides an independent cross-check of our picture. 

The ratio of the energy density in the shock-injected 
vorticity (or shear flow) to the thermal energy density is 
given by 



Eshcar ^.TTLpV^ 



3fcT2 



0.14, 



(13) 



where kT2 — 2.4 keV and uj_ ~ 400kms~ . Interest- 
ingly, this value is very close to the turbulent-to-thermal 
pressure support of £tuib/£th< — 0.2 at i?200 found in cos- 
mological simulations of the formation of galaxy cluste r 
(|Rvu et al.ll2008l: iLau et"aI|[200llB"attaglia et al.ll2010D . 
We caution that our estimates strongly depend on our 
model assumptions for the orientation and position of 
the shock surface with respect to the LOS since projec- 
tion effects of a circular torus could account for a fraction 
of the observed ellipticity, reduce the shear contribution, 
and increase the inferred curvature radius. Another com- 
plication is the amount of magnetic helicity inside the 
torus that acts as a stabilizing agent against the shear 
that t ries to distort and potentially disrupts the plasma 
torus (jRuszkowski et al.ll2007t IBraithwaitell2010D . 

5. SCRUTINISING OUR MODEL OF NGC 1265 

5.1. Examining the Physics of the Model in Detail 

Using results on the shock properties and geometric pa- 
rameters from Sects. [31 and m we demonstrate the phys- 
ical feasibility of this model. 

5.1.1. The radio bubble and torus 

The expansion time of the bubble in the fllament (prior 
to shock passage) is of order the external sound crossing 

^ The sign of the post-shock vorticity is given by the barochnic 
term, Vp X VP, which in the top view of Figure 1 (lower panel) 
points out of the page. 
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thue, Texp — R/ci ~ 0.4 Gyr, usmg a a sound speed 
of the pre-shock region ci ~ 285 km s~^ (fcT/0.4 keV)!/^ 
and a bubble radius of i? ~ 150 kpc. This assumes 
that the outer radius R is unchanged upon transfor- 
mation into a torus (as suggested by numerical simu- 
lations of En filin fc Brii ggcn 2002) and that R equals 
the minor axis of the ellipse of the projected torus that 
is unaffected by the shear. We can estimate the ener- 
getics of the previous outflow by calculating the PdV 
work done in inflating the bubble. Using the pre-shock 
gas pressure of Pi ~ 3.6 x 10"^'' erg cm~'^, we obtain 
i^bubbio = l/(7roi-l)Pil^ ~ 4x 1058 erg using 7,.^! = 4/3. 
This fits well inside the range of observed jet energies of 
FRl sources and implies an equipartition magnetic field 
in the bubble of Bi,eq = (87rPi/2)i/2 ~ 0.66 //G that 
becomes after shock passage -82, oq = Si^eqC"^/^ ~ 3 /zG 
with C ~ 10; somewhat higher than i?eq — 1 ob- 
taine d by equipartition argume nts applied to the radio 
flux (jSiibring fc de Bruvnlll998l ). The radio synchrotron 
radiating electrons of Lorentz factor F emit at a fre- 
quency i/gyn given by Equation ([Ij. Using a range of 
magnetic field values of _Boq — (1~3) /iG, we require elec- 
trons with a post-shock Lorentz factor r2 ~ 4800 — 8500 
to explain the observed radio torus at 600 MHz and de- 
rive a minimum value of Fi = Fs/C^/^ ~ 2200 - 4000 
for the upper cutoff of the electron population in order to 
be able to account for the observed radio emission after 
shock passage. Synchrotron and inverse Compton aging 
of fossil electrons in the bubble occurs on a timescale 



67rTOeC 



' syn, ic, 1 



fTT(S; 



mb 



^i,oq)ri 



= (0.6 - 1) Gyr > t,-_ 



(14) 

Calculating the maximum cooling time that the bubble 
could have been hibernating as given by Equation ()14p . 
we show that we meet this criterion for Fi so that the 
model bv lEnfihn fc Gopal-Krishnal ()200lD applies. 

5.1.2. The current outburst and radio tail of NGC 1265 

Assuming that the twin tails together constitute a 
cylinder of radius rtaii — 1.'5 and projected length 
Lt, tail — 11', we can estimate the volume for the ra- 
dio tail of NGC 1265, Vtaii ^ Trr^^iiLt.taii/tane* ~ 3.8 x 
10™ cm"^ usin g 9 :^ 32°. Taking the minimum energy 
density from Sij bring fc de Bruvnl (jl9 98) in the tail of 
NGC 1265 to be around £eq — 10~^^ erg cm"'^ we esti- 
mate the energy content of the radio tail to be about 
Ejot — 3.8 X 10^8 erg which nicely coincides with the 
result that we derived for the previous outburst. The 
necessary accreted mass to power the outflow amounts 
to Mace — Ejct/ivc^) — 2 X lO^Mo, where we adopted 
a typical values of 77 ~ 0.1 for the efficiency parame- 
ter. Using the time since shock crossing in our model 
of T ~ 1.8 X lO^yr, we estimate the jet power of about 
6.7 X 10^^ erg s"^ which is in range of observed jet lu- 
minosities of FRI sources. Using the 1.4 GHz fiux for 
NGC 1265 of 8 Jy, we obtain a radio luminosity of 
vL^ ~ 7.5 X 10^° erg s~^, or a 1% radiative efficiency 
which is plausible. 

5.1.3. Estimates from the jet bending 

Using the high-resolution V LA observations of the 
head structure of NGC 1265, lO'Dea fc OwctI ([l986l ) 



give a minimum jet pressure /jet,min — (1 ^ 3) x 
10~^^ ergcm"'^ which is about 10 times the ICM esti- 
mate at the present position of the galaxy, PicM(''gai) — 
1.4 X lO^^^ergcm"^ with rgai ~ 1.45 Mpc (see Sects. [3] 
andHJ. The jet radius is rjet 1" (360 pc), whereas 
the projected bending radius of the jet is around 30" (11 
kpc), so the physical ratio is rh/rjot ~ 80 when account- 
ing for deprojection and finite resolution effects. Using 
Equation (jA2p . this leads to a Mach number ratio of 
7Mjot/-Mgai — 2.8. The temperature estimate at the po- 
sition of the galaxy is fcT(rgai) ~ 3.1 keV, so that we 
obtain A^gai — 2.8 and eventually A^jot — 7.8. If we take 



the jet power to be ijet/2 Pjetv 



(where Ljot is 

the luminosity of the 2 jets and we assume that most of 
the jet power is kinetic energy), the power can be ex- 
pressed in terms of the jet pressure and Mach number 
as 

r 

7jot-PjotA^j^ct^^jot7rr2^f (15) 



^JCt 



Since our estimate of Lj^t — 6.7 x lO'*^ erg s ^, we get 
Vjct — 7000 km s~ which is a typical jet velocity. Solv- 
ing for the number density of the jet material, we get 
njot — Ljct/(27rTOpW?,^r?,j) ~ 10~'^cm~'^ which is a fac- 
tor of three larger than the surrounding ICM, ?^(7'gal) ~ 
3 X 10~^cm~'^ and again a plausible value. 

5.2. Stability Analysis of the Torus 

There are three mechanisms that have the potential to 
destroy the torus once it has formed; namely, (1) large 
scale shear flows injected at the shock front, (2) ICM 
turbulence, and (3) Kelvin-Helmholtz instabilities at the 
interface due to the vortex flow around the minor cir- 
cle of the torus. By considering each of these processes 
separately, we will show that none of them can seriously 
impact the torus on an eddy turnover timescale of the 
stabilizing vortex flow, 



'^cddy 



^2,11 



2 X 10® yr. 



(16) 



where rmin — 25 kpc denotes the minor circle of the torus, 
^^2,11 — V cos cj)/Cs — 740 km s^^ is the post-shock gas ve- 
locity parallel to the shock normal, v ~ 2550 km s~^ the 
total velocity of the galaxy, and ~ 9° is the shock 
obliquity defined in Table [T] Here and in the following, 
we apply the numerical values of our consistent model 
as derived in Sects. [3]and|4l On a timescale larger than 
Teddy, dissipative effects will thermalize the kinetic en- 
ergy of this stabilizing flow, and eventually external shear 
flows and turbulence might start to destroy the torus de- 
pendi ng on the magnetic helicity and morpholog y in the 
torus (jRuszkowski et al.ll2007t lBraithwaitell2010l) . 

(1) The amount of vorticity Wghear injected at a curved 
shock with curvature radius i?curv is given by (jLighthilll 
|1953) 



l^she 



(C, -1)2 vx_ [Cs-lfvsim 



Cs 



Rc 



Rc 



(17) 



where Cs is the shock compression factor and v± is the 
(post-)shock gas velocity perpendicular to the shock nor- 
mal. Note that for vorticity injection into an irrotational 
flow encountering a shock, one necessarily needs a curved 
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shock surface according to Crocco's theorem ()1937[ ). The 
generation of vorticity rehes on a gradient of either the 
entropy or stagnation enthalpy along the shock front 
(perpendicular to the shock normal). This can be easi- 
est seen by considering a homogeneous flow encountering 
the shock surface at some constant angle (obliquity). It 
experiences the same "shock deflection" and amount of 
entropy injection along the shock front. This is because 
mass and momentum conservation across the shock en- 
sures the conservation of v± and yields z;2,|| = Vi \^/Cs- 
If the bubble was impacting the shock surface at larger 
angle, the induced stabilizing vortex flow around the mi- 
nor circle of radius r^^in is accordingly smaller since the 
perpendicular velocity component cannot contribute to 
the vortex flow, 



"2,11 



(18) 



We can then derive a criterion for the stability of the 
torus in the presence of shock injected shear flows, 
i^shcar < i^torus that Can bc cast into a requirement for 
the shock obliquity. 



< 'pent = arctan 



Rn 



.(a - 1) 



80°. 



(19) 



Here Cs — 3.4, i?curv — 850 kpc as inferred from the 
observed ellipticity of the torus, and r^^i^ ~ 25 kpc. Since 
our preferred value of the shock obliquity cj) — 9° 
0crit , the torus can easily be stabilized by the vortex flow 
around its minor circle against large-scale shear. 

(2) To assess the impact of ICM turbulence on the 
torus' stability, we compare the turbulent energy density 
Eturb to the kinetic energy density of the torus vortex 

flow £vort, 



Eturb 



Evort 



Eturb 



£th.2 



Eturb / Ri 



Eth,2 
Evort 



Eth,2 



-2/3 



3fcT2 



0.04.(20) 



Here, eth.2 and kT2 — 2.4 keV denote the post-shock en- 
ergy density and internal energy, fi = 0.588 denotes the 
mean molecular weight for a medium of primordial ele- 
ment abundance, and nip denotes the proton rest mass. 
To quantify the ratio Eturb/Evort in our model, we as- 
sume a turbulent pressure support relative to the ther- 
mal pressure at i?200 of around 20% when correct ing for 
the energy contributi o n due to bulk flows (.Rvu et al.l 
12001 ILau et al.l &M IBattadia etalWwuh . We fur- 
thermore assume that the turbulent energy density is 
dominated by the injection scale (which should be com- 
parable to the curvature radius at the shock). Using 
a Kolmogorov scaling of the turbulence, the power per 
logarithmic interval in scale length is down by another 
factor of (i?curv/''min)^^'^ — 10 at the minor radius rmin 
compared to the injection scale i?curv Hence, the en- 
ergy density of the stabilizing vorticity flow around the 
minor axis of the torus outweighs any turbulent energy 
density of the ICM on this scale by a factor of about 
25! Equation d^U]) enables us to solve for the maximum 
shock obliquity allowed so that the torus is stable against 



turbulent shear flows. 



< arccos 



Eturb / R: 



Eth 



-2/3 



3fcT2C| 



1/2 



78° 



(21) 

where v ~ 2550 km s ^ and Cs — 3.4. Interestingly, 
this is a similar value as i^crit obtained in Equation (jl9p 
and is another manifestation of our flnding that Eturb ^ 
Eshcar (Equation (|13p ) implying that a large fraction of 
post-shock turbulence could have been injected at curved 
shocks. 

(3) We flnally show, that the Kelvin-Helmholtz insta- 
bility at the interface of the torus due to the vortex 
flow should be suppressed by internal magnetic fields. 
The shock compression acts mostly perpendicular to the 
forming toroidal ring. Consequently, the field component 
parallel to the bubble surface and hence to the surface of 
the torus is preferentially amplified and dominates even- 
tually the magnetic energy density (Enfilin & Briigge^ 
[200l . Within the torus, we assume equipartition mag- 
netic fields of i3eq,2 — (l-3)/zG which correspond to ra- 
dio and pressure equipartition values, respectively (Sec- 
tion [53J.® Assuming the torus to be in pressure equi- 
librium with the surroundings and a relativistic electron 
component in the bubble whose pressure is dominated 
by momenta ~ mgC, we can derive a density contrast 
between the radio plasma and post-shock density of the 
ICM, 

"torus _ _ ^ ,, ,„_3 



5 X 10" 



?^2,ICM 



where kT2 — 2.4 keV. We can now compare the ratio of 
the Alfven velocity, ua, in the surface layer of the torus 
to the velocity in the vortex flow around the minor circle 
of the torus, V2^\\ , 



VA 

"2,11 



Br,, 



11, 



AirSp. 



(23) 



2,ICM^2,|| 



where n2,icM = 2 x 10^"^ cm^'^. In such a configuration, 
magnetic tension is strong enough to sup press the Kelvin- 
Helmholtz instab ility of the vortex flow (jChandrasekhad 
[T9MtlD^[2n07h . To summarize, if the magnetic field is 
not too far from equipartition values, B > i?2.cq/10, then 
the line of arguments presented here guarantees stability 
of the torus for at least an eddy turnover timescale Teddy 

6. DISCUSSION AND CONCLUSIONS 

In providing a quantitative 3D model for NGC 1265, 
we provide conclusive observational evidence for an ac- 
cretion shock onto a galaxy cluster. The accretion shock 
is characterized by an intermediate-strength shock with 
a Mach number of ~ 4.2l°'2. In the context of the 
ICM, it appears that even weak to intermediate strength 



° If the magnetic field was sub-equipartition in some regions 
just below the torus surface, the external vortex fiow implies a 
matching internal vortex flow that acts in piling up the field in a 
layer just internal to the surface. In principle internal vortex flows 
could be on small scales as long as the effective internal vorticity 
matches the external one. However buoyancy of the magnetic field 
in combination with an effective reshuffling of the magnetic fleld 
due to internal vortex flows should give rise to a surface layer flUing 
magnetic fleld. 
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shocks of ^ 2.3 are able t o accelerate electrons a s 
the radio relic in A521 suggests ()Giacintucci et aLll2008[ ). 
This supports the view that these intermediate strength 
shocks are very important in understanding non-thermal 
processes in clusters; in particular as they are much 
more abundant than high-Mach num ber shocks (e.g., 
iRvu et all 120031 iPfrommer et all 120061 ) . The presented 
technique offers a novel way to find large-scale forma- 
tion shock waves that would likely be missed if we only 
saw giant radio relics. These are thought to be due 
to shock-accelerated rel ativistic electrons at merging or 
accre tion shock waves (jEnfilin et al.l 119981 : iBagchi et al.l 
|2006| ). whereas the discrimination of the two types of 
shock waves is difficult on theoretical as well as observa- 
tional grounds (|Enfilinii2006) . 

The rich morphology of the large scale radio structure 
that arches around steep-spectrum tail of NGC 1265 en- 
ables a qualitative argument for the need of shear to 
explain the orientation of its ellipsoidally shaped torus. 
The energy density of the shear flow corresponds to a 
turbulent-to-thermal energy density of 14% — consistent 
with estimates in cosmological simulations. However, 
simulations are needed to quantify the effect of magnetic 
helicity internal to the bubble on the torus morphology. 
In any case, the presence of post-shock shear implies 
the amplification of weaker seed magnetic fields through 
shearing motions as well as the generation of magnetic 
fields through the Biermann battery mechanism if the 
vorticity has been generated at the curved shock. Even 
more interesting, the shock passage of multiphase gas 
with a large density contrast necessarily yields a multi- 
ply curved shock surface that causes vorticity injection 
into later accreted gas on the corresponding curvature 



scales. Hence, shearing motions are superposed on var- 
ious scales which implies amplification of the magnetic 
fields on these scales. 

More sensitive polarized and total intensity observa- 
tions of NGC 1265 are needed to confirm our model 
predictions with radial polarization vectors of the radio 
torus and a polarization degree of around 5%. In addi- 
tion, future observations of different systems similar to 
NGC 1265 would be beneficial to get a statistical mea- 
surement of properties of the accretions shocks and its 
pre-shock conditions. This provides indirect evidence for 
the existence of the warm-hot IGM. Our work shows the 
potential of these kind of serendipitous events as ways 
to explore the outer fringes of clusters and dynamical 
features of accretion shocks that are complementary to 
X-ray observations. 
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APPENDIX 



A. DERIVATION OF THE JET CURVATURE 

Although fairly standard, for completeness we show the derivation of the jet curvature radius that is caused by an 
external ram pressure wind due to the motion of the galaxy through the ICM. We denote the mass density, velocity, 
and radius of the jet by pjet, fjot, and rjet, respectively. The two jets coming out of the active galactic core back 
to back are assumed to be initially a cylinder of length /jet- Each jet is then bent over a bending radius rf, by the 
ram pressure wind of mass density and velocity picM and v. We equate the jet momentum /9jotWjot7rr?,jZjct with the 
transverse force due to the ram pressure wind that acts over a jet propagation timescale (along the bended path in 
steady state), picMW^2rjot/jctTb/(2ujot), to obtain the following equality 

PlCMV = Pjct^^jct^-jct TT- (Al) 

ZV^ct Z 

Solving for the ratio of bending-to-jet radius, we obtain 

''jet -^gal 7ICM^'iCm' 

where we introduced the Mach numbers of the jet and the galaxy, TWjot — Wjot/cjot and A^gai — f/cicM, the adiabatic 
exponents of jet and surrounding ICM, 7jot = 4/3 and 7icm = 5/3, and their pressures, Pjet and Picm- 

B. DERIVATION OF THE TRANSFORMATION TIME OF A SHOCKED BUBBLE 

When a low-density bubble crosses a shock of speed Wsi in the ICM, the shock accelerates into the bubble, pulling 
post-shock ambient gas with it. The original bubble-ICM CD follows the shock at a speed intermediate to that of 
the incident and bubble shocks when the incident shock is at least moderately strong. Because the shock intrusion 
begins sooner and impacts more strongly on the leading edge of a round bubble than its periphery, the ambient gas 
penetrates the center of the bubble first. An initially spheroidal bubble will then evolve into a torus (vortex ring) on 
a timescale determined by the crossing time of the CD through the bubble. 
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Fig. 4. — Left: ratio of the speed of a shock inside a low d ensit y bubble, Dgb, to the incident shock speed, Vgi, for density contrasts, 
5 = 10~^, 10~^. Solid curves represent solutions to Equation IIB5I1 . Dotted curves represent the internal bubble sound speed in units of 
the external shock speed in the I CM . Right: ratio of the contact discontinuity speed inside the bubble to the incident shock speed from the 
relevant expression in Equation ||B2|| . 



The shock-produced dynamics in the bubble can be estimated from the exact solution to the simple ID Riemann 
problem of a plane shock impacting a CD. This is found most simply in the initial rest frame of the bubble, which we 
assume to be in pressure equilibrium with the unshocked ICM and has a density, pb = Spi, with S < 1. At impact a 
forward shock will penetrate into the low density bubble with speed, Vsb, while a rarefaction will propagate backward 
into the post-shock ICM. The bubblc-ICM CD will move forward at the same speed as the post-shock flow in the 
bubble. The full Riemann solution is obtained by matching the pressure behind the forward shock inside the bubble 
to the pressure at the foot of the rarefaction in the ICM. 

We define the ICM and bubble adiabatic indices as 7i and 7b respectively. Assuming the shock is propagating from 
the left, we have right to left four uniform states, 0, 1, 2, 3, where and 1 are separated by the forward shock, 1 and 2 
are separated by the CD, while 2 and 3 are separated by the reverse rarefaction (note that this numbering scheme differs 
from the main body of the paper). The state "0" represents the initial conditions in the bubble, while "3" represents 
conditions in the ICM post-shock flow. We set Pi ~ P2 = P* and vi ^ V2 ~ v^. — ucd- The initial bubble (and ICM) 
pressure is Pq, while the initial bubble and ICM sound speeds are related by ci, — cq ~ \/^bPo/po = Q \/"/b/iliS)- The 
Mach numbers of the external, ICM shock and the internal, bubble shock are Aii — Vsjci and M.b = Vsh /cb = p^Aii 
and we introduced the Mach number ratio p.^ 



From standard shock jump conditions, we have 



^2 7i - 1 



^3 = ^(7.-M^-^)Po 



(7. + l)M^ ^ 1 (7. + l)M^ 
(7,-l)X? + 2^' <5(7.-l)X? + 2' 



--4t^- (B1) 

,2a.2 76-1 



7fc + 1 V 27 
_Jjb + lWMl_ 
{^b^l)p^Mj + 2P' 

2 [p'M^, - 1) 

-^CD =Vi=V2^V^ = — — Cb. 

-fb + 1 pMi 

In addition, the Riemann invariant connecting states 2 and 3 through the rarefaction, along with the relation p cx P'^/T' 

® Note that the variable has a different meaning compared to the main part of the paper where it denotes the mean molecular weight. 
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inside the rarefaction give 



while 



= V3 + 



-C3 



1-1^ 
P3 



7.P3 _ 2 [t.A^I - V] [(7. - l)Ml + 2] 



These can be combined to give the following equation for fi — Mb/ Mi- 



(B2) 



(B3) 



1 



7. + 1 M^Xf - 1 

76 + lM^l/2(^2_i) 



1/2 



1 [2^,Ml-i. + lf^[{i..-l)MWA 



,1/2 



li 



1 



Mi 



1 



(7» + l)(27bM'>''-7b + l) 

(7fc + l)(27,;A^2_^^^l) 



(B4) 



which can be solved numerically for /i, given Mi, 7i, 76, and 5. For i5 
our problem we expect 7^ = 5/3 and 7^ = 4/3. This gives 



1, 7i = 76 the obvious solution is /x = 1. For 



1 = 3 U 



1/2 



m'X? - 1 



[(5Mf-l)(Xf + 3)] 
1 



1/2 



V 5A^?-1 



1/5' 



(B5) 



Equation (|B5p is displayed for (5 = 10^^ and 5 = in the left panel of Figure |H 

In the strong (external) shock limit, ^iMi 3> 1, it is easy to see that approximately /i oc \/6. Empirically we find in 
this limit when 5 that /x ^ 2\/5, so that Vsh ~ 2usi. As Mi — > 1, the solution /i — )■ 1, so Mb 1 also applies. The 
more general solutions for small to moderate incident shock strengths are shown for two values of 6 as solid curves 
in Figure m The lower bound for Vsh is given by Cf, = Vsi \/ 76 j^i [l/(A^i-\/5)] and is shown in each case by a dotted 
curve. We note that the solution shown in Figure |4] is almost identical to the case of 7^ = 7b, since the internal shock 
in the bubble is barely supersonic with A^f, 1 due to the large sound speed in the bubble. This can be easily seen 
by taking the ratio of the corresponding solid-to-dotted lines which provides Mb for each external Mach number Mi- 

We can understand the general behavior of an increasing shock speed inside a low-density bubble relative to the 
incident shock speed for smaller Mach numbers Mi by the following line of arguments. For small Mach numbers, 
both waves are just nonlinear sound waves, so each propagates near the local sound speed, which is much larger in 
the bubble. It turns out for large Mach numbers Mi that the pressure just behind the penetrating shock is smaller 
than the external post-shock pressure by a factor that scales with the density contrast b. This pressure drop comes 
from the rarefaction going back into the shocked ICM. For weak shocks that rarefaction is weak (provided the original 
bubble was in pressure equilibrium). 

The right panel in Figure S] shows the behavior of wcd corresponding to the shock solutions in the left figure panel. 
In the limit /x^A^f 1, the expression for tico in Equation (|B2I) takes the form I'cD/'fsi ~ (2/ (7b + l))\/7b/^7i 
Note that the limit n^Mi ^ 1 requires the internal bubble (pscudo-) temperature not to be too high and hence 
implies a lower limit on the bubble density contrast of (5 > 10^"^ assuming the original bubble was in pressure 
equilibrium. Applying our empirical result for strong shocks with 10"'^ < (5 <C 1 (namely, /i ^ 2\/J), we would 

expect vcd/vsi — >■ (4/(7b + l))v7f)/7i ~ 1-5, which is within about 20% of the exact solutions in this limit. This value 
for vcd/vsi is also a reasonable estimate at moderate shock numbers, say Mi > 3. For Mi > 2 the limit ?;cd > fsi 
still applies. As Mi approaches one, however, the speed of the CD drops below the incident shock speed, since as the 
bubble shock becomes a sound wave, the CD's motion vanishes. The results of these considerations are quantitatively 
confirmed by a suite of 2D axisymmetric simulations where we varied the Mach number Mi — {2, 3, 5, 10} and initial 
bubble/ICM density contrast (see Figure [2] for one realization). 

For our moderate-strength shock example in Perseus we can use these results to estimate the time for the CD to 
cross a bubble diameter, 2i?bubbic, transforming the bubble into a vortex ring, with a toroidal topology; namely. 



2i?b 



ubblc 



1.4 X 10** yr. 



(B6) 



Here we adopted our fiducial values of i?bubbic — 150 kpc, (5 ~ 5 x 10 ^ (Equation ([22]) ). and Wgi = MiCi — 

~ 1400 km s"^ where we used the sound speed in the pre-shock gas in 
0.4 keV and fi^^ = 0.588. We emphasize that this estimate is close to the 



Mi^/^kTi/{fj,mwmp) ~ 4.2 x 330 km s 
the infalling filament of temperature kTi 

value r3^4 = 1.2 x 10^ yr that we adopted in our model and consistent with the derived the error bars on Mi 
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